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Abstract 

A fractional Hamiltonian formalism is introduced for the recent com- 
bined fractional calculus of variations. The Hamilton-Jacobi partial differ- 
ential equation is generalized to be applicable for systems containing com- 
bined Caputo fractional derivatives. The obtained results provide tools to 
carry out the quantization of nonconservative problems through combined 
fractional canonical equations of Hamilton type. 

MSC 2010: Primary 26A33, 49K05; Secondary 49S05 
Key Words and Phrases: fractional canonical formalism, Hamiltonian 
approach, variational principles of physics, nonconservative systems, com- 
bined fractional derivatives, variational calculus 



1. Introduction 

The conservative physical systems imply frictionless motion and are 
simplification of the real dynamical world. The systems in Nature always 
contain internal damping and are subject to some external forces that do 
not store energy and which are not equivalent to the gradient of a poten- 
tial. For nonconservative dynamical systems the energy conservation law 
is broken and, as a consequence, the standard Hamiltonian formalism is 
no longer valid for describing the behavior of the system. Moreover, not 
only the energy, but also other physical quantities such as linear or angular 
momentums, are not conserved: Noether's classical theorem ceases to be 
valid pi]. 

In order to model nonconservative dynamical systems, a novel approach 
to variational calculus was proposed by Riewe in 1996 and 1997, [2l"l [25] . 
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The idea is based on the concept of fractional differentiation and frac- 
tional integration, [H[9]. Roughly speaking, the appearance of fractional 
equations is natural due to physical reasons of long-range dissipation, frac- 
tional wave propagation in complex media and long memory. Fractional- 
order Euler-Lagrange equations and the Hamilton formalism with frac- 
tional derivatives are investigated in |16[ [23] . For a fractional Hamiltonian 
analysis of irregular systems, see |5]. 

The fractional calculus of variations is nowadays an emerging field - 
see [U El [7J Uni HZ] and references therein. A general fractional variational 
calculus, based on a combined Caputo derivative, was recently introduced 
by the authors in [19] and further developed in [20] and [21]. In [3], a 
generalization of the fractional variational calculus is given by considering 
generalized Hilfer's fractional derivatives with two and three parameters. 
Here we study the importance of the new combined variational calculus in 
fractional mechanics, describing nonconservative dynamic systems. 

In Section [2] we recall the combined fractional Euler-Lagrange equa- 
tions. Our results are then given in Section [3l In Section \3. II we show that 
combined fractional Hamiltonian equations of motion can be obtained for 
nonconservative systems. The constants of motion cease to be valid and a 
new notion is introduced in Section [3.21 Canonical fractional transforma- 
tions of first and second kind are studied in Section 13.31 and Section 13.41 
respectively. Subsequently, in Section 13.51 an Hamilton-Jacobi type equa- 
tion is derived and a fractional quantum wave equation suggested. 

The understanding of the Hamiltonian dynamical systems, classical 
or quantum, has been a long standing theoretical question since Dirac's 
quantization of the classical electromagnetic field. The combined fractional 
Hamiltonian approach here introduced seems to be a promising direction 
of research. 



2. Preliminaries 



Let a, j3 £ (0,1) and 7 6 [0,1]. The fractional derivative operator 
c Dj'^ was introduced in Malinowska and Torres [19], as 



7^ + (l- 7 )^\ 



Cn/3 



7 



For 7 = and 7 



erf 



Dl and °Z> 



C — C na 



1 we obtain the standard Caputo operators: c Dq" 



aD". Consider a fractional Lagrangian 



L[t,q, C D^ 



I) 



depending on time t, the coordinates q = [q%,. . . , Qat] and their fractional 



velocities c D^' p q 



71 



91, 



c n»N,0i^ 



The Euler-Lagrange 
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equations corresponding to (pQ) were studied in [20] and form the following 
system of N fractional differential equations: 

^ + D^ ^^ = 0, i = l,...N, (2) 

dqi d c D^Qi 

where 

with a Dx' and X D^ denoting the classical left and right Riemann-Liouville 
fractional derivatives. 



Remark 2.1. With the notation from [2"0"1 EZTj. 

dL dL 
di+iL = — and d N+1+i L = , i = l,...N. 

dqi d c D^' IJl q i 

In contrast with [20^ 121] . we adopt here the notation used in mechanics 

[hies]. 



More general necessary optimality conditions than ([2]) are found in 
|20}l21j. However, all our previous works [19 |, I20 |. [2Tj on combined fractional 
derivatives consider the Lagrangian approach only. In contrast, here we 
develop the combined fractional Hamiltonian formalism. We show that ([2]) 
can be reduced to a special Hamiltonian system of 2A fractional differential 
equations: the combined canonical fractional equations. 



3. Main Results 



In analogy with the classical mechanics, let us introduce the canonical 
momenta pi by 

dL 

Pi = a—, i = l,...N. (3) 

3.1. Combined Euler Lagrange equations in Hamiltonian form 



Assume that 



d(pi,...,p N ) 



'7i 91.— > u in In j 

function theorem, we can locally solve equations ([3]) with respect to 



7^ 0. Then, by the implicit 



- L/ 7i yi j 



c Dj£'^ N qN. The fractional Hamiltonian is defined by 



N 



(4) 
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where ^-D^ are regarded as functions of variables t, q±, . . . , </at>Pi, 
Pn- Therefore, the differential of H is given by 

V dH , A dH , dH 



From the defining equality we can also write 

N N N 

dH = Y, C D^q t dtp l + Y J V l <fB^q l - £ — dq t 
i=l i=l i=l ™ 



V — d c D ai>/ V - — dt (6) 



C D^q t , ^- = D^tPi, i = l,...,N, (7) 



The terms containing dF D^l q, t in ^ cancel because of the definition of 
canonical momenta. Applying relations ([2]), we get 

n N fit 

dH = ^ C D^ gi dp t + £ D^Pidqi - °-gdt. 

i=l i=l 

Comparison with ([5]) furnishes the following set of 2N relations: 

dpi ll *' dqi 

which we can call the combined fractional canonical equations of Hamilton. 
They constitute a set of 2N fractional order equations of motion replacing 
the Euler-Lagrange equations (|2]). Moreover, 

dH dL 

For integer-order derivatives, it can be shown that 

dH _dH _ dL 
~dt ~ ~dt ~ ~dt 

(see, e.g., [TH p. 220]). Hence, if L (and, in consequence of (jHJ), also H) 
is not an explicit function of t, i.e., in the autonomous case, then H is a 
constant of motion. This is a well-known consequence of Noether's theorem 
[271 [29] • With noninteger-order derivatives this is not the case, see [T2| [13] . 
Observe that using equations ([7]) , we can write 

dH _ v ( c n a *>*a dpi + n^ n dqi \ + dH 

1=1 v ' 

In contrast with the integer-order case, the terms containing fractional and 
classical derivatives of coordinates and momenta do not cancel. Therefore, 
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in general we obtain nonconservative systems and classical constants of 
motion cease to be valid 124] . To deal with the problem we introduce 
the notion of fractional constant of motion. 

3.2. Fractional constants of motion 

To account for the presence of dissipative terms, we propose the follow- 
ing definition of fractional constant of motion. 



Definition 3.1. We say that a function C(i,q, p) is a fractional 
constant of motion of order (a,/3,y) if D^f 1 [t \- > C (i, q(i), p(t))] is the 
null function along any pair (q, p) satisfying the 2N combined fractional 
canonical equations (|7|). 

It follows from ([7]) that if qi is absent in the fractional Hamiltonian, then 
DflZ l Pi = 0, i.e., if = 0, then pi = w . — is a fractional constant of 

motion of order (a^, 7$). 

3.3. Canonical fractional transformations of first kind 

We now look for transformations under which the combined fractional 
canonical equations of Hamilton ([7]) preserve their canonical form. Consider 
the simultaneous transformation of independent coordinates and momenta, 
qi and Pi, to a new set Qi and Pj with transformation equations 

Qi = Qi(*,q,p), Pi = Pi(t,q,p), i = l,...,N. 

The new Qi and Pi are canonical coordinates provided there exists some 
function K(t,Q,P) such that the equations of motion in the new set are 
in Hamiltonian form: 

dPi ~ ^ Qu dQi ~ u ~ ' • • • ' [y) 

Transformations for which equations Q are valid are said to be canonical. 
Function K plays the role of the Hamiltonian in the new coordinates set. 
If Qi and Pi are to be canonical coordinates, then they must satisfy the 
modified fractional Hamiltonian principle of form 

5 [ (E Pf^Qi ~ K(t, Q, P)^ dt = 0. (10) 
At the same time, the original coordinates satisfy the similar principle 

5 f a ^p/^'Ni-tf^p)! dt = 0. (11) 
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The simultaneous validity of equations (jlOp and (jlip means that integrands 
of the two integrals can differ at most by a total derivative of an arbitrary 
function F, sometimes called a gauge term [28] : 

, N N 

jF + Y, PfD*f l Q* ~ * (*> Q- p ) = E« C£, £ - H & * p)- 

8=1 1=1 

The arbitrary function i 7 works as a generating function of the transforma- 
tion, and is very useful in developing direct methods to solve problems of 
the calculus of variations [18] and optimal control [261 EO] • Function F is 
only specified up to an additive constant. In order to produce transforma- 
tions between the two sets of canonical variables, F must be a function of 
both old and new variables. For mechanics involving fractional derivatives 
(cf. [M]) we need to introduce variables % and Qi, i = 1, . . . , N, satisfying 



dt ~ ^ QU dt ~ ^ 



For 7j = 1, aj = 1, = 0, i = 1,...,N, these new coordinates are 
the same as the usual canonical coordinates. However, when dealing with 
fractional derivatives, the coordinates % and Qi will not be canonical, so all 
canonical expressions must be written in terms of c D".'^qi and c Dy.'^*Qi. 
For a generating function F\ (t, q, Q) , the integrands of (fTOl) and (fTTI) are 
connected by the relation 

, N N 

-Fx (t, q, Q) = E Pi C D^ qi - H(t, q,p)-£ PfD^&Qi + tf(t, Q, P). 



eft 

Because 
we obtain 



i=l i=l 



,r, v^^Fi , ^ n dF x , 8Ft , 

£=1 



(12) 



" ^< » = -fls' - Pi= 9Q"' (13) 

i = 1, . . . , N. We call to (|13p the canonical fractional transformations of 
first kind. 



3.4. Canonical fractional transformations of second kind 

We shall now introduce canonical transformations of second type with 
a generating function Fz(t, q, P). We begin by noting that the transfor- 
mation from qi and Qi to q~i and Pi can be accomplished by a Legendre 
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transformation. Indeed, by equations (fl~3j) . —Pi = J^. This suggests that 
the generating function F 2 can be defined by 



N 



P 2 (t, q, P) = Fx (t, q,Q) +^ P;Qi- 



From equation (|12p we obtain 

i=l i=l 

, / AT 

' , F 2 (t,q,P)-J]P f Q i 



i=i 



i=l i=l 



Hence, 



-H)dt+[ Y,Qi^T + K)dt = dF 2 (t,q,P). 



dt I \ ^— ' 

\i=l / \i=l / 

Repeating the procedure followed for Pi in Section 13.31 we obtain the 
transformation equations 

-gf = K-H, Pl = —, Ql = —, » = !,...,*. (14) 

3.5. Fractional Hamilton Jacobi equation 

Under the assumption that K is identically zero, for integer order 
derivatives we know, from equations ([U]), that the new coordinates are 
constant. The same situation occurs, under appropriate assumptions, in 
the fractional setting when in presence of Riemann-Liouville or Caputo 
derivatives |15[ [22] . In our general combined setting this is not the case, in 
particular for the new momentum P. To solve the problem we proceed as 
follows. Assume that the transformed Hamiltonian, K, is identically zero. 
Then the equations of motion are 

11 = ^-^ = 0, H = ^;;^ = o, i = i,...,N. 



The function K is related with the old Hamiltonian H and the generating 

dF 

dt 



function P by K = H+%£. Hence, K will be zero if, and only if, P satisfies 
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the equation H(t, q, p) + ^ = 0. It is convenient to take F as a function 
of cji and Pj. Then, from equations (|14p . we can write 



/ dF 2 \ dF 2 



Equation (I15p is a fractional version of the Hamilton-Jacobi equation. 
Therefore, the quantum wave equation with the combined C D^'^ fractional 
derivative is suggested to be 



In this paper we give an Hamiltonian formulation to the general frac- 
tional Euler-Lagrange equations obtained in |19j . Motivated by the recent 
results in the literature of Physics dealing with fractional mechanics and 
nonconservative systems [TJ [151 122] , we illustrate some possible applications 
of our combined variational calculus in mechanics and quantization. 
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